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ABSTRACT 

The X-ray spectra of accreting black hole systems generally contain components 
(sometimes dominating the total emission) which are well-fit by thermal Comptonization 
models with temperatures ~ 100 keV. We demonstrate why, over many orders of 
magnitude in heating rate and seed photon supply, hot plasmas radiate primarily by 
inverse Compton scattering, and find equilibrium temperatures within a factor of a few 
of 100 keV. We also determine quantitatively the (wide) bounds on heating rate and 
seed photon supply for which this statement is true. 

Plasmas in thermal balance in this regime obey two simple scaling laws: 
Ott — 0.1(lh/l3y^^', and a ~ 1.6(Zs/Zh)^''^. Here the hot plasma heating rate 
compactness is Ih, the seed photon compactness is la, the temperature in electron 
rest mass units is G, and the Thomson optical depth is tt- The coefficient in the 
first expression is weakly-dependent on plasma geometry; the second expression is 
independent of geometry. Only when h/lh is a few tenths or greater is there a weak 
secondary dependence in both relations on tt- 

Because a is almost independent of everything but h/lh, the observed power law 
index may be used to estimate h/lh- In both AGN and stellar black holes, the mean 
value estimated this way is h/lh ~ 0.1, although there is much greater sample dispersion 
among stellar black holes than among AGN. This inference favors models in which the 
intrinsic (as opposed to reprocessed) luminosity in soft photons entering the hot plasma 
is small, or in which the hard X-ray production is comparatively distant from the source 
of soft photons. In addition, it predicts that Qtt ~ 0.1 - 0.2, depending primarily 
on plasma geometry. It is possible to construct coronal models {i.e. models in which 
Is/lh — 0.5) which fit the observed spectra, but they are tightly constrained: tt must be 
~ 0.08 and 9 ~ 0.8. 



1. Introduction 

Ever since the paper by Shapiro et al. (1976), it has been a commonplace to explain the hard 
X-ray spectra of accreting black hole systems, whether they are galactic binaries or AGN, by thermal 
Comptonization. Indeed, such models often give very good descriptions of the broad band X-ray 
spectra of these systems (Sunyaev & Titarchuk 1980; Grebenev et al. 1994; Johnson et al. 1994; 
Zdziarski et al. 1994; Haardt et al. 1993). In recent years ART/P, SIGMA and OSSE observations 
have substantially increased the volume of data available for constraining these models, but the basic 
framework has been sustained: for all systems identified on the basis of independent arguments as 
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accreting black holes, thermal Comptonization is an acceptable description of the X-ray spectrum 
from ~ 3 - 300 keV to within the uncertainties in both measurements and theory. 

The best fit values for the model's two free parameters, temperature and optical depth, are 
almost always in the range 30 - 300 keV and 0.1 - 5, respectively (Haardt et al. 1993; Bouchet 
et al. 1993; Johnson et al. 1993; Zdziarski et al. 1994). There is certainly no selection effect that 
could explain this comparatively narrow range of parameters in the AGN, for most are discovered by 
optical, ultraviolet, or radio searches. While galactic black hole candidates are almost all discovered 
by X-ray techniques, these searches would have been sensitive to a much wider range of temperatures: 
for fixed luminosity and optical depth, many of these sources could easily have been detected if their 
temperature had been as low as 1 keV or as high as ~ 1 MeV. 

In most previous work in this field. X-ray spectra fit with thermal Comptonization models 
are described in terms of the inferred temperature and optical depth. For many purposes this is a 
reasonable procedure, for the spectral formation models are most directly framed in terms of these 
parameters, and they do have a direct physical meaning. However, if we are to understand why these 
conditions are found, we will do better to speak in terms of heating rates, cooling rates, and dynamics. 
In lower temperature plasmas, studies of how dissipative heating and photoionization balance 
recombination, bremsstrahlung, and collisional excitation led to the realization that the dissipation 
rate and ionization parameter determine the temperature and ionization state. In high temperature 
plasmas, dissipation and Compton recoil or pair production using externally-created high energy 
photons balance annihilation, bremsstrahlung, and inverse Compton scattering of low energy photons; 
it is the primary goal of this paper to learn how this balance scales with the externally-imposed 
parameters and boundary conditions. 

A principal result of this program will be to explain why it is that inverse Comptonization 
dominates other radiation processes in these systems, and why its parameters are restricted to such 
a comparatively narrow range. We will argue that if one frames the issue in terms of the physically 
causal parameters, namely the heating rate and the seed photon supply rate, the observed parameters 
are approximate "fixed points" of the nonlinear equations which determine them. This approach 
also turns out to have a second important consequence: when attention is restricted to values of 
temperature and optical depth which are actual equilibria, simple scaling laws emerge relating 
observables to the causal parameters. 

We begin by outlining the physical processes which operate in these systems (§2), and then 
describe how we compute equilibria (§3). In §4 we show quantitatively over how wide a range of causal 
parameters inverse Comptonization dominates the cooling rate. We also show that the equilibrium 
temperature does not depend in any direct way on the traditional Compton j/-parameter, but that a 
different scaling relation does describe the equilibrium quite well. In addition we discover a simple 
relation between the power-law index and the ratio of soft compactness to hard. We then discuss in 
§5 observable consequences of the scaling relations, and why the causal parameters for accreting black 
holes almost always lie within the (large) volume of parameter space corresponding to the narrow 
range of observed parameters. We summarize in §6. 



2. Thermodynamics of Hot Plasmas 

To highlight the basic physics of these hot plasmas, we focus on a highly idealized model 
(following, e.g. Bjornsson & Svensson 1991, BS91 hereafter): a spherical, homogeneous cloud of pure 
ionized hydrogen in which all particle distribution functions are (relativistic) Maxwellians at the 
same temperature O = kT/rne(? (see Pietrini & Krolik 1994, PK94 hereafter). This temperature is 
maintained by a heating rate per unit volume H = 3Lh/{^nR^), where R is the radius of the cloud 



-3- 



and Lh is its luminosity. Photons with initial energy xome(? {xq <C O) are supplied externally in such 
a way that the photon intensity — at xq and all other photon energies xrUeC^ — is isotropic and constant 
at every point in the cloud. 

Two sorts of equilibria must be established: between pair production and annihilation; and 
between heating and cooling. Because the cooling rate depends on the number density of leptons, 
while the rate of pair production depends on the number density of high energy photons, these two 
equilibria are strongly coupled. Although (approximate) fully self-consistent solutions have been 
obtained for the version of this problem without an independent source of soft photons (BS91), 
previous treatments of equilibria in which externally generated soft photons enter the Comptonizing 
cloud (Zdziarski 1985; Haardt & Maraschi 1991, 1993; Ghisellini & Haardt 1994) have a number of 
limitations, some of which we remove in the present work. For example, Haardt & Maraschi (1991, 
1993) have modeled the production of X-rays in AGNs with a system in which a hot corona of total 
{i.e. electron plus positron) Thomson optical depth tt < 1 covers an accretion disk. In their model 
the entire accretion power of the disk is dissipated in the corona. Because the disk thermalizes and 
reradiates the majority of the X-rays directed toward it, the ratio between the seed photon luminosity 
Ls and the corona luminosity Lh is fixed at approximately 0.5. We extend the range of possibilities 
to include both larger optical depths and other ratios between the seed photon luminosity and the 
plasma luminosity. In addition, their thermal balance does not include annihilation losses, internal 
photon production in the corona, or the Comptonization losses on the internally-generated photons. 
In their limited range of parameters, these approximations are justified, but we wish to study these 
physical circumstances over a wider range of parameters. On the other hand, Ghisellini & Haardt 
(1994) cover a much larger range in the luminosity ratio Lg/L^, but they assume a pure pair plasma, 
and so are limited to follow a narrow track in heating rate and optical depth. 



2.1. Key parameters 

Lightman (1982) and Svensson (1982) showed that it is a very good approximation to suppose 
that pair balance equilibria depend on the net lepton optical depth Tp rather than on the net lepton 
density or size of the plasma cloud separately. Note that in the pure hydrogen limit, Tp = UparR, 
where rip is the proton density. 

Because Svensson (1984: S84) also showed that the number of pairs z = n^/Up, where n+ 
is the positron number density, could be found as a function of G and Tp, most treatments since 
have taken as one of the independent variables of the system. In fact, for a given fixed value of 
G, the pair balance equation, h+ = fiA (where n+ and tla are the total pair production rate and 
the pair annihilation rates, respectively), can be solved to give z — z(Q,Tp), independently of the 
thermal balance of the system, i. e. , without taking care of how much heating is actually required to 
maintain this temperature. Moreover, with this approach, one finds that there are either two distinct 
equilibrium solutions or no solutions at all for a given choice of (B, Tp), thus implying the existence of 
two equilibrium "branches" for the equilibrium curve z = z{Q) for a given Tp value (see S84, BS91, 
PK94). One branch corresponds to solutions for which z 1, and is determined essentially by photon 
(pair-creation) processes; the other branch is controlled (mostly) by particle processes, and z ^ 1. In 
a, z = z{Q) plot the two branches connect at a maximum allowed value for the plasma temperature, 
Qmax{Tp); at greater temperatures the pair production rate always exceeds the annihilation rate, so 
that no equilibrium can be found. 

Once the pair balance has been found, for a given 0, the heating rate H (per unit volume) 
required to maintain that temperature can be determined by equating it to the sum of all the cooling 
rates, which can be written as functions of G, Tp, and z. In causal terms, however, the temperature 
is a derived parameter which depends on H, therefore we can choose to take H as the parameter 
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and derive consequently the temperature, together with z, by coupUng the thermal balance ccjuation 
to the pair balance one. As pointed out in BS91, the thermal balance equation can be written most 
concisely in terms of the compactness 

= (1) 
me.(? 3 

Thus, for a given value of r^, this equation and the pair balance equation may be inverted to find 
0('hiTp) and z{lh,Tp). These solutions are single-valued, and, therefore more useful than the ones in 
terms of ©. 

Our work differs from that of Svensson in adding a new element, external soft photon input. 
Two more parameters (or more) are required to specify this. We suppose, for simplicity, that the 
soft photon spectrum is a 5- function at energy xo- We call their (constant) number density no; the 
corresponding dimensionless quantity is the ratio no /tip. It is useful to also treat the density of soft 
photons in terms of the alternative dimensionless quantity, the soft photon compactness, 

^^ = ^^' (2) 

meC^ R 

where Ls is the luminosity of soft photons necessary to maintain the density no- The actual value of 
Ls depends on geometry; we assume that it is the luminosity which exactly balances the rate at which 
soft photons leave the cloud: 

r° 

Ls=4nR' diJ.iJ.Is=4.n'^R'^Is, (3) 



where the intensity of soft photons is given by 



laix) = -^meC^xnoS{x — Xo). (4) 
47r 

Thus, the soft photon compactness definition (eq. (2)) reduces to 

Is = TTTpXo . (5) 

Up 

In summary, the equilibrium conditions of our system are determined by four key causal 
parameters: 

1) cloud (hard) compactness lu; 

2) cloud net lepton (or, equivalently, proton) Thomson optical depth r^; 

3) external soft photon compactness h (or no /rip); 

4) soft photon energy xo- 

For a given choice of these, from the solution of the coupled pair and thermal balance equations 
one derives the other key parameters: 

a) plasma temperature = Q{lh,Tp,ls,xo), 

b) pair density, z = n+ jnp =S> total optical depth to scattering tt = aT{'2n+ + np)R = {2z + l)rp; 

We stress that over the very large volume of parameter space we consider solutions parameterized 

by Is/lh do not depend strongly on xq. This is because Comptonization is unsaturated (in the sense 
that over the great majority of the parameter space of interest, most of the injected photons are not 
scattered all the way up to the electron temperature), so that the Compton cooling rate depends more 
directly on the energy density of seed photons than on their number density. Therefore, we can focus 
our attention on the first three key parameters mentioned above and their interrelations. 
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2.2. The Compton Parameter y 



It has long been conventional to describe the results of Comptonization by a single parameter, 
which we call 

y„r = rciax{TT , tt) , (6) 

where maj(;(TT,Tr) is an estimate of the mean number of scatterings a photon suffers before escape, 

and AO is the the mean amplification factor, that is the average fractional photon energy increase 
due to a scattering event with an electron drawn from a Maxwellian distribution at non-relativistic 
temperature G (when the electron energy is much larger than that of the photon). This ynr 
parameter (the subscript nr means non relativistic) arises as the key dinicnsionless quantity when one 
writes a Fokker-Planck equation for the evolution of photon energies as a result of Comptonization 
(Kompaneets 1957; Rybicki & Lightman 1979). 

In certain commonly encountered limits, it has several important physical interpretations. In the 
limit in which <C 1 and tt ^ 1, the shape of the emergent spectrum may be written in terms of 
y„r (lUarionov & Sunyaev 1975, Sunyaev & Titarchuk 1985). When <C 1 and y„r <Si ln(3©/xo), the 
average ratio between the photons' energy on escape and their energy upon entry is simply exp{ynr) 
(lUarionov & Suuyaev 1975; Shapiro et al. 1976). Sometimes {e.g. Liang 1979) this is taken to mean 
that the ratio of the emergent luminosity to the seed luminosity is exp(j/„,.), but this is not correct. 
As pointed out by Sunyaev & Titarchuk (1985) and Dermer et al. (1991), the output luminosity is 
usually dominated by photons on the tail of the energy change distribution, not the photons whose 
energy change is near the mean (see also §4.4.2). It also takes more time than average for these higher 
energy photons to emerge because they must scatter more times to reach those energies. In addition, 
when the output spectrum is "hard" enough, the photon energy range which dominates the luminosity 
approaches the energy G, so that saturation effects can be important for them, even if they axe not 
important for the average photon. For all these reasons, exp{ynr) is not a very good expression for the 
relation between Lh and Lg. In §4.4.3 we derive a much more accurate scaling law. 

There are other difficulties which also limit the utility of the y-parameter for the problem we 
consider. When is not small compared to unity, photon energies do not change by small fractional 
amounts per scatter, and a Fokker-Planck description is therefore of questionable validity. As for the 
relation of y with the spectral shape, this turns out to be quite problematic as well; in fact, even when 
tt ^ 3, we often deal with transrelativistic or mildly relativistic temperatures, bringing back the issue 
just discussed above. On the other hand, we are also interested in cases in which tt is not very large, 
and for values of tt ~ 1 or smaller, the diffusion approximation breaks down, and the spectral shape 
is not directly related to y (see § 3.1.1), even in principle. 

Despite these difficulties, several attempts have been made to find generalizations of the 
t/-parameter to relativistic temperatures. All the different attempts agree about how to generalize the 
expression for the mean photon amplification per scattering to relativistic temperatures: 

,1 _i I ^3(1/9) _ Ki(i/e) , 

A._l + 40^^^^^Q^ -l + 40^^^^^^^ +160 , (7) 
where Kn is the nth-order modified Bessel function, and, following S84, its usual approximation is 



=i 1 + 40 -h 160^, (8) 



at least for not too large values of 0. 



However, there is significant disagreement about the best way to generalize the y parameter. 
Loeb, McKee, & Lahav (1991) make the hypothesis that the probability distribution for the number 
of scatterings that a photon can suffer is a Poisson distribution. They then find that the ratio of 
the average energy of the escaping photon to its initial value is, similarly to the non-relativistic case, 
exp(j/tot), where the parameter ytot = {Ar — 1) < u >, and < w > is the average number of scatterings 
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for a photon before escape. Therefore, in this view, the representative relativistic Compton parameter 
can be defined as yr = ytot, so that, approximating Ar with eq. (8), 

yr ~ (49 + 16e^)TT (l + y) . (9) 
where < u >~ rT(l + tt/3) is also a reasonable approximation. 

On the other hand, Zdziarski et al. (1990) estimate the ratio of escape energy to initial energy 
is by supposing that every photon escapes after suffering the mean number of scatters. With this 
approximation, the relativistic Compton parameter is 

J/. =ln(l + 4e + 16e^)rT(l + rT/3). (10) 

where the degree of approximation for Ar is the same as in eq. (9) . 

In the range of physical conditions for the Comptonizing plasma wc arc interested in, the two 
different generalizations discussed above give rather different estimates for the parameter, and, in 
practice we find that neither of them really directly helps describing even qualitative features of 
Thermal Comptonization in the mildly relativistic temperature regime, as we discuss in § 4.4.2. 

As a matter of fact, in this regime, for oven more reasons than those discussed above for 
the general CEise, the use of a Compton parameter for a description of the conditions of thermal 
Comptonization is not much justified. This is probably true for what regards the characterization of 
the spectral shape as well, for which the relevance of the Compton parameter is actually still an open 
question (see Titarchuk 1994, and Zdziarski, private communication). 



3. Calculation 



As we remarked before, we must solve two equilibrium conditions, one for pair balance and one 
for thermal balance. Here we describe the physical processes we included, the approximations we 
made, and how we arrived at our solutions. 



3.1. Pair balance 



We write the pair balance equation as follows: 

riA = (n+)77 + {n+)ie + {n+)ip, (11) 

where the annihilation rate per unit volume tia is equated to the sum of the relevant pair production 
rates, namely photon-photon (indicated by the subscript 7 — 7), photon-electron (7 — e), and 
photon-proton (7 — p) . Particle-particle pair creation rates are neglected here, since they are always 
negligible (see PK94 and references therein). We also assume that no pairs escape from the cloud. 

All of the pair creation rates included depend on the total Comptonized radiation field in the 
plasma cloud, which is given by the sum of the contribution due to Comptonization of the internally 
generated photons, and that coming from Comptonization of the externally created soft photons. For 
this reason, we must define a representation for both contributions. 
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3.1.1. Comptonized Radiation Field 



Following S84, we approximate the contribution from Comptonization of the photons originating 
in the cloud plasma with the photon number density spectral distribution 



nc{x) = nw{x) + npix) 



2 02 e 





■ e-^/e" 


+ 


ni 




X 



(12) 



nc{x) is the sum of a Wien spectral component, nw{x), and a "flat" component, npix). The factor 
ni is the flat component normalization, explicitly defined in BS91; n-y is the number density of Wien 
photons from Comptonization of internally generated radiation, and it is derived as 

n-f = tescifsriBr + focriDc), (13) 

where fiBr and hoc are the energy integrated photon production rates for bremsstrahlung and for 
Double Compton process respectively, explicitly defined in Svensson (1984) and in PK94; /b and foe 
( again defined in S84) are the probabilities that a soft photon of bremsstrahlung or Double Compton 
origin is scattered into the Wien peak before escaping the cloud; these also depend on the parameter 
yi defined in S84 which describes the approach to Wien equilibrium for internally-generated photons. 
Therefore, ni and n-y are functions of Q and tt- The photon escape time from the cloud, tesc, is a 
useful quantity to give a rough description of radiative transfer effects, and it is here defined as in 
BS91, 

te.e(rT,e) = (I) (l+iflrT-T) , (14) 

where gT{Q)TT = Tw is the scattering optical depth including the Klein-Nishina cross section decline 
at high energies, by averaging its effects on a Wien photon distribution, and gT{&) is (from BS91) 

^^^"^ ~ \ (3/16)e-2(ln(1.12e) -|-3/4)[l + (0.1/9)]-^ for 9 > 1 ' 

the factor 1/3 is related to the spherical geometry chosen (see BS91). 

There is another radiation field related relevant quantity that must be determined, and it is 
the energy Xm below which photons get strongly absorbed; this, in general, represents the lower 
limit for spectral integrations. As a matter of fact, this energy is defined through a relation which is 
coupled to the pair and thermal balance (since it depends on 9 and z) and must be actually solved 
simultaneously with those conditions (see PK94). 

Now we consider the contribution to the radiation field inside the cloud due to Comptonization 
of the soft photons originally created by an external source. Since we are analyzing the problem in 
the framework of the thermal pair equilibrium model originally studied by S84 and BS91, we can be 
satisfied with a description of the Comptonized external soft photon contributions to the radiation 
field which is at the same level of approximation as the one we use for the Comptonized internally 
generated photons. A representation of this type is given by Zdziarski (1985, herafter Z85; 1986): the 
Comptonized "soft photon" spectrum is given, in terms of spectral photon number density nsc(x), as 
the superposition of a power law with an exponential cutoff at the plasma temperature 9 and a Wien 
peak at the same temperature, 

n,e(x) = no,|-(Qj —- + -^—^[^-) —-^=n,,{x)+ns„(x). (15) 

In this relation uq^ represents the Wien photon density from Comptonization of external soft photons, 
and nop is a normalization factor for the power-law component. The relation that best defines 
the spectral index a of the power-law component in terms of 9 and tt depends on tt (Z85, and 
corrections in Zdziarski 1986): 



-8- 



For tt larger than a few, a is well described by a rolativistic generalization (Z85) of the 
non-relativistic form given by Sunyaev & Titarchuk (1980): 

1/2 



a= (4+7(TT,e)j --, 
where 7(rT, 6) is given by Z85, 



7(tt, 6) 



3(Tr + 2/3)2effe ' 
with the relativistic correction 

5e = Kz{l/Q)/K2{l/@), 

again given in Z85. 

• For small tt, the best description is a generalization by Z85 to general Q and tt 1 of the 
formula originally derived by Pozdnyakov et al. (1977) for cases with tt <^ 1 and 3> 1: 

a = -^, (16) 

where 

^ = l + 4ese, (17) 

and 

PrT=i-^ [2r^ - 1 + e-'"" (27T + 1)], (18) 
(corrected for a misprint, as explained in Z86). 

As for the ratio (now/nop), this is given by Z85 as 

now _ r(a) 



nop r(2a + 3) 
which is again a generalization of ST80 results. 



Pr^irr), (19) 



We have now to relate the normalization factor nop to the total number density of incoming 
external soft photons, no. From the conservation of photon number, in the absence of absorption 
processes, we have that 



nod - e--) = ^\.(.)d. = ^ {/r + (^) /2} ^ 



where 



(20) 



Ixa/e 

and I2 — 2. The lower limit of integration chosen, xq/Q, implies that we neglect recoil effects in 
the Comptonization process. The integral Ji can be calculated in terms of Whittaker functions, 
following Gradhsteyn & Ryzhik (1980) (integral 3.381.6), and then, since Xo/Q << 1, analytically 
approximated, by using expressions from Abramowitz & Stegun (1972) (par. 13.1.33, 13.5.5-13.5.12), 
so that we finally obtain 

which is valid for any value of a (< or > 1), and 

\np J \np J h +2{now/nop) 

When <C 1 and tt S> 1 (where the ST80 approximations are valid), these approximations are 
fairly good fits to the ST80 solutions, particularly at the low and high energy extremes. However, they 
underestimate the ST80 photon intensity by factors of a few in the neighborhood of a; ~ 0. When 
TT <C 1, these approximations agree reasonably well with Monte Carlo simulations (Zdziarski 1986). 
Pair production rates computed on the basis of these approximations should always be fairly accurate 
because they are most sensitive to the highest energy photons. 



3.1.2. Radiative Pair Production Rates 



Having made these choices for the representation of the contribution to the radiation field in the 
cloud due to Comptonized external soft photons, we now describe how we calculated the corresponding 
contributions to photon-photon and photon-particle pair production rates. 

All pair production rates involving photons make use of the Wien component, so our general 

procedure is to substitute + no^ for in the appropriate expressions in Appendex B of S84. This 
is a simple substitution for the photon-particle pair production rates, and for all those contributions 
to the photon-photon pair creation rate except the interaction of Wien photons with the power-law 
part of the external photon Comptonized spectrum. For this last part, from the general integral form 
given by S84, we use the expression 

2 r'^ 
{h+)u,pi = ^nopK + non,)e"-^ J <P{so)s-^"+^^^^K^+3{2sl^ye)dso, (22) 

where 4>{so) is a function given by Gould & Schreder (1967), with corrections from Brown, Mikaelian, 
&: Gould (1973). Thus, the pair creation rate depends on the strength of Comptonization, so 
that the pair balance, thermal balance, and emergent spectrum calculations must be performed 
self-consistently. 



3.2. Thermal Balance 



In a schematic form, the equation of heat balance is 

^Ai = H, (23) 

i 

where each separate cooling mechanism is denoted by Aj. The processes we include are: 

a) Bremsstrahlung cooling, Ab, which is the sum of the rates for e — p, e — e, and e"'"e~ interactions; 

explicit expressions are given in BS91. Here we must mention that there is one more mechanism 
of internal photon production, the double Compton process (7-l-e— >7-|-7 + e). This can be 
important in the "low" -temperature range (6 < 0.1), in terms of number density of generated 
photons (S84, PK94). On the other hand, from the energetics point of view, it is negligible 
with respect to other cooling mechanisms, since the created photons are generally at a much 
lower energy than the original interacting photon, usually one of the Comptonized high-energy 
photons (see PK94, Thorne 1981) 

b) Annihilation losses, defined hero as Aa ~ An — 2mec'^nA, whore Aa is the total annihilation 

cooling rate; this definition of the cooling rate accounts for the fact that we are interested in 
the losses to the thermal energy of the plasma, and the rest mass energy loss contribution is 
actually balanced by pair creation, since pair density is in equilibrium. 

c) Inverse Compton losses on all (i.e. bremsstrahlung and, possibly. Double Compton) internally 

created photons, written as 

(Ac)6r- = meC 3Qn~f/tesc, 

where n,, is given by eq. (13). 

d) Finally, Inverse Comptonization of the external soft photons, including formation of both the 

Wien component and the power law component: 

1) The Compton cooling term due to the contribution of external soft photons to the Wien 
peak is 

(Ac)sio = meC^30noi„/tesc, (24) 
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2) The Compton cooling due to the power law portion can be evaluated as 

2 /• oo poo 
m-eC / f \ J _ 2 nop ^ a -ti. _ 2 "Op 



(Ac)p! = — — / npi{x)xdx = rrieC ^ 8 / t "e dt = nieC ^ 8Jo(a, a:o/9), 



xo/e 



(25) 



where 



„oc r r(l - a) for < a < 1 

Io{a,xo/e)= r"e"'dt ~ <^ (x„/e)-i/2[|ln(a;o/e)| - V(l)] for a = 1 .(26) 

■^^o/e l(x-o/e)('-"V(a-l) forl<a<2 

In this expression ^{z) is the logarithmic derivative of the gamma function, and 
ip{l) = 0.577156649. 

The escape time tesc is given by eq. (14). 

One final correction is necessary: the preceding expressions actually give the total emergent 
spectrum, which sums the inverse Compton cooling rate with the luminosity injected in soft 
photons. The net Compton cooling due to processing of the external seed photons is therefore 

(A,), = (A,)3» + iAc)pi - ^J^i^ - exp{-TT)]. (27) 



4. Results 



4.1. Independence with respect to xq 



In the following subsections all the cases we illustrate were computed with the same initial soft 
photon energy, namely xo = 10~^. To check the sensitivity of our results to the choice of xo, we have 
also computed equilibria with xq a factor of 10 larger and smaller. On the basis of these results, we 
can assert that, provided the basic condition O/xo 3> 1 is fulfilled, for a fixed value of Is/lh there is 
only a very weak dependence of the equilibrium solutions on the actual value of xo- For an order of 
magnitude change in xq, the corresponding variations in the equilibrium values of O and tt typically 
amount to only a few per cent and are never larger than ^ 15%. We are therefore confident that our 
results are, indeed, representative of the general problem. 



4.2. Pair equilibrium curves: the effects of changing no/up 



Wc begin by presenting a global view of our equilibria, with an emphasis on illustrating how 
a variable soft photon supply can affect them. In order to facilitate comparison with previous work 
{e.g. S84, BS91, PK94), in which G was generally taken as the independent variable. Fig. 1 displays 
the normalized pair density z = n+frip, the total optical depth tt, the required heating rate (as 
a compactness), and the spectral index a of the output Comptonized power law component of the 
spectrum, all as functions of G. These equilibria share the same Tp (= 1) and used our standard value 
of xo, 10~^. In each panel we show a family of curves, parameterized by no/up. We remind the reader 
that fixed no /rip is equivalent (for constant TpXo) to fixed Is. 

Increasing no /rip clearly has a number of dramatic effects on the character of the equilibrium: 
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a) The maximum temperature, Qmax, that the system can attain steadily decreases; indeed, in 

almost all cases, the temperature that can be achieved for a given heating rate Ih decreases as 
Is increases. 

b) At fixed temperature, the pair content increases on the \ow-z branch, and decreases on the high- 2: 

branch. 

c) As corollaries of the previous point, at fixed temperature on the high-2 branch tt decreases and 

a increases, while the behavior is precisely opposite on the low-2 branch. 

d) Again for fixed G, on the low- 2: branch increasing no/rip implies an increase in 1^; however, when 

z > 1, Ih becomes almost independent of no/ Up for fixed and the different curves Zft(0) of 
Fig. Ic all converge. 

e) At very large z (and therefore Ih), the z{Q) curves become independent of h- The value of Ih at 

which this happens increases linearly with no/ Up; that is, the curves converge at a single value 

of la/lh- 

We now briefly elaborate on these points in order; most will be re-discussed from the point of view of 
Ih as the independent parameter in the following subsection, and the physics behind them will then 
become much clearer. 

That the temperature falls as no /np increases is no surprise; the number of photons available for 
Compton cooling increases, forcing the temperature to decrease. Note that the range of temperatures 
permitted also narrows as no/up increases. 

Similarly, when there arc many soft photons available, fewer electrons arc required to account 
for the total cooling at fixed temperature, so the equilibrium pair content falls on the high-z branch. 
On the other hand, there is a larger number of photons which may potentially be scattered to energies 
above the pair production threshold, so the number of pairs on the low- 2: branch increases. 

Fewer pairs translates directly into smaller tt, and smaller tt implies less Comptonization. 
Therefore, increasing no/ Up spells a softer power-law on the high-a branch, and a (very slightly) 
harder one on the low- 2 branch. 

The increase in Ih at fixed G on the low- 2; branch noted in point d) is a corollary of point a). 
Similar reasoning explains a related effect, not illustrated in Fig. 1: increasing Tp decreases Qmax 
at fixed no/ Up. With a fixed number of seed photons, but greater opportunity to scatter them, the 
temperature falls. For example, with no/rip = 10^, Qmax falls from 1.12 for Tp = 0.1 to 0.216 for 

Tp = 1. 

The insensitivity of the lh{Q) curves to Is on the high-2 branch is due to a combination of 
different effects. For fixed Is, the limit of very large Ih produces a condition in which 2 S> 1, so that 
Tt becomes very large and Comptonization is saturated. The output spectrum is then dominated, 
of course, by the Wien portion. In that limit, the ratio between Wien photons and particles in the 
plasma is a function only of temperature because pair balance determines the number density of 
photons. When this is the case Ih also depends only on G (S84). 

At values of Ih slightly lower than those at which complete Wien equilibrium is attained (i.e., for 
slightly higher temperatures on the high-2 branch), at fixed G, larger Is reduces 2 and hence tt (on 
the low- 2; branch, Is has little influence on tt, of course). This in turn leads to a reduction in the 
amount of energy carried off by the average escaping photon because Comptonization is weakened. 
However, this reduction is almost exactly compensated by the increase in the number of photons 
available due to the increase in Is. We will be able to describe this phenomenon more clearly and 
quantitatively armed with the formalism and scaling laws developed in §4.4. 

However, as Fig. 1 shows, when Is > 3, this compensation is no longer effective. In this case, 
the compactness curve crosses the bundle of curves describing smaller Is equilibria, and does not join 
them until true Wien equilibrium is reached at substantially larger values of Ih (see also Fig. 1 in 
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Svensson 1986 for another ilhistration of this effect). When Is is this large, the greater number of 
photons available for Comptoriization becomes more important than the reduction in the mean energy 
per photon, and a larger heating rate is required to maintain thermal balance. Put another way, more 
energy is required to approach Wien equilibrium because there are more photons to scatter up to 
energies ~ kT. 

Point e) is due to the fact that when z ^ 1, the increase in total density leads to a 
large increase in the internal photon production rate by brenisstrahlurig. At sufficiently large z, 
internal photon production dominates over external soft photon input and the Compton cooling on 
these internally generated photons becomes, in turn, comparable to or dominant over other cooling 
contributions. At this point the hard compactness is basically determined by the "internal" conditions, 
i.e., it is (almost) independent of the external photon input. In Fig. 2 we show an example of this 
behavior, for the case {Tp = 1, no /rip = 10"*). The various cooling rates Aj are plotted in the form of 
compactnesses: 

_ 47r A, 

3 ^"nieC^npic/R)' ^ 
where the subscript i indicates the specific process. As Fig. 2 shows, over most of the range of Ih, 
inverse Compton cooling on externally-created photons accounts for nearly all the cooling, but at the 
highest values of Ih, it is overtaken by inverse Compton cooling on internal photons. An inspection 
of the cooling rates for other cases confirms that the onset of this regime (and the consequent 
convergence of the z{Q) curves) corresponds to the condition in which the compactness ratio h/lh 
reaches a specific (rather extreme) value, namely ~ lO""*. 



4.3. Plasma Temperature as a Function of Ih, Is/lht and 



As we have just seen, discussing our equilibria in the conventional language using Q as the 
independent parameter makes many effects difficult to explain because they require arguments based 
on self-consistency that appear circular. Much greater clarity is achieved by transforming to Ih as 
independent variable because fundamentally, the temperature is the result of equilibration between a 
given heating rate and the various cooling mechanisms. Similarly, while tt is more closely connected 
to the degree of Comptonization than Tp, the extra number of electrons and positrons which changes Tp 
to Tt is a consequence of the pair equilibrium, while Tp is (we presume) controlled by the forces acting 
on the plasma. Finally, for the remainder of this paper we will use h/lh as the third independent 
variable because it is more directly related to the global energy budget (see §4.4) and the geometry of 
the plasma than is no /rip: 

where Ls^intr is the intrinsic luminosity of the source of soft photons, C is the fraction of solid 
angle around the hot plasma occupied by optically thick matter capable of absorbing the X-rays 
and re-emitting their energy in lower energy photons, d is the typical distance from the soft photon 
source to the hot plasma, and 4> accounts for possible transfer effects within the hot plasma. Because 
Is/lh has such a simple physical interpretation, it has been widely utilized elsewhere {e.g. Haardt & 
Maraschi 1991, 1993; Ghisellini & Haardt 1994; Dermer et al. 1991). 

In this section, we present the results of our calculations in terms of these three independent 
parameters, and describe the trends we find. In the following section we will explain these trends 
physically. The results we show here span a large volume of parameter space, but with the greatest 
dynamic range in lu: 0.03 < < 3 x 10^; 0.01 < Is/lh < 1; and 0.1 < Tp < 3. (We have actually 
computed equilibria for Is/lh = 0.001 as well, and the corresponding results are included in the figures 
referred to in § 4.4.3 and in the related discussion.) This range of parameters should cover most cases 
relevant to the study of accreting black holes (see §5 for more discussion of this point). We chose 
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{Tp)max = 3, because, with increasing Is, higher values of Tp would yield temperatures which fall below 
our range of interest (for example, for Tp = 3 and Is ^ 0.94, Qmax — 0.11). 

In Figs. 3,4,5 the equilibrium curves as functions of Ih are shown for h/lh = 1.0,0.1,0.01 
respectively; in each figure we plot 0{lh), together with the curves for ttHh), ct{lh), and yr{lh) (as 
defined in eq. (9)). Each quantity is also computed for several values of Tp. In Fig. 3 {h/lh = 1), 
Tp = 0.1, 0.5, 1.0, 2.0; in Figs. 4 and 5, Tp = 0.1, 0.5, 1.0, 3.0. 

The most prominent feature of these plots is the wide region in Ih over which, taking Tp and 
Is/lh fixed, Q is very nearly constant. When Is/lh ^ 0.01 and Tp > 0.1, the temperature changes by 
no more than a few percent over as much as four decades in heating rate! 

The existence of this constant temperature regime is not new; it is implicitly contained in the 

results of Haardt & Maraschi (1991). It appears (sec §4.4 for a fuller discussion) whenever inverse 
Compton scattering of externally supplied photons is the dominant cooling mechanism, Is/lh is fixed 
(as, for example, by absorption and subsequent reradiation in soft photons of a fixed fraction of the 
hot plasma's radiative output), and there arc few pairs, so that tt — Tp. The new points we make here 
are to define the range of conditions in which it applies, and to elucidate the physics which controls it. 

The actual value of G on the flat part of the curve is also rather insensitive to h/lh'. an order 
of magnitude increase in that quantity decreases the temperature by only a factor of 2 - 2.5. It is, 
however, more sensitive to changes in Tp: in the range of values we have explored, d\nQ/d\nTp ~ — 1 
if Ih is flxed somewhere in the range for which the temperature is constant (Haardt & Maraschi 1993 
find a similar dependence of G on Tp for the specific case of h/lk — 0.5). 

The dynamic range in Ih corresponding to constant temperature also depends on both the 
compactness ratio h/lh and Tp. For fixed Is/lh, larger Tp yields a wider range of Ih in which the 
temperature is constant. Likewise, for flxed Tp, increasing Is/lh also leads to a widening of the interval 
in hard compactness over which G is constant. 

The upper boundary in Ih for the constant temperature range can be anywhere from ~ 1 to 
> 3 X 10^. It is set by the point at which pair density starts to be no longer negligible {z > 0.05), 

so that Tt begins to become greater than Tp, and this is sensitive to both h/lh and Tp. In those 
cases where G remains flat up to the largest value of Ih shown, that is because z remains very small 
throughout that range of Ih ■ 

The position of the lower boundary varies rather less than the position of the upper boundary. 
Throughout our parameter range it is ~ 0.01 - 0.1. As we shall show in detail in the following 
subsection, the lower boundary is set by the point at which other cooling processes (notably 
bremsstrahlung) become competitive with inverse Compton cooling. The position of this boundary 
depends hardly at all on Is/lh, but does depend somewhat on Tpi increasing Tp increases Ih at the 
lower boundary. 

At the very highest values of Ih, all the equilibrium curves corresponding to different values of 
Tp, for a given ratio h/lh, tend to coincide, again going towards Wien equilibrium conditions for the 
plasma. This is clearly seen in Fig. 5, corresponding to the smallest Is/lh shown here. This is in fact 
another manifestation of the independence of G with respect to Ih when tt is fixed; it is just that 
when Ih is very large, z is also very large, and tt now depends on Ih rather than on Tp. In fact, if we 
had plotted these results for fixed tt, the constant temperature range would extend to even greater Ih, 
and the different curves would remain separate even at very large Ih. Note that the larger the value 
of Is/lh, the larger the value of Ih at which this occurs. This is because larger Is/h leads to lower 
temperature, and consequently, for fixed Ih, a lower rate of pair production. 

We close this section by pointing out that, for fixed Is/lh, ot is independent of Ih over an even 
broader range of Ih than G is; in addition, it is rather less sensitive to tt. Several factors contribute 
to this. In the range of Ih for which G is constant, tt is also constant, so a, which is a function of 
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those two variables (§3.1.1), must certainly be constant. However, at larger Ih, where tt grows with 
the greater importance of pairs, O falls with increasing If^ in such a way as to cancel the effects of 
changing tt (see §4.4.3). Similarly, the weak dependence on Tp in the \ow-z regime is due to the 
inverse relation of O and Tp (also elaborated more fully in §4.4.3), while the weak dependence on Tp in 
the high-« regime is due to the decoupling of tt from Tp when z ^ 1. The not result is that, provided 
Ih is greater than the lower boundary of the constant temperature range, a is very nearly a function 
of le/lh alone. 



4.4. Cooling Physics 



4.4- i- Relative importance of specific mechanisms 



To explain the trends described in the previous subsection, we now examine the relative 
importance of the different cooling mechanisms. To illustrate this, we plot the ratios of the different 
cooling compactnesses, h, to the total hard compactness (i.e., the heating rate compactness) Ih as 
functions of Ih itself for the same values of h/lh and Tp shown in Figs. 3, 4, and 5. Thus, Figs. 6, 7, 
and 8 show these plots for Is/lh = 1.0, 0.1, 0.01 respectively, with each of the four panels corresponding 
to a given value of Tp. The curves shown in each plot are: the ratio Ihc/lh of the total Compton cooling 
compactness to the total compactness; the ratio Ib/Iji of the bremsstrahlung cooling compactness to 
the total compactness; and the ratio Ia/Ih, where I a is the annihilation cooling compactness. 

It is apparent that Comptonization losses, Ihc, including both the contribution due to 
inverse Compton on the internal bremsstrahlung photons and that of inverse Compton on the 
externally-supplied soft photons (that is, 

^--f;;^:^3^[(A^V + (A.).l, (30) 

with the notations of section 3.2), completely dominate over a large range in total hard compactness. 
Compton cooling dominates all the way from the range where lu is so large that pair processes 
dominate down to the bottom of the constant temperature range. 

Ultimately, at sufficiently small Ih, bremsstrahlung overcomes inverse Compton cooling. The 

value of Ih whore this occurs (and therefore the lower boundary to the constant temperature range), 
lh,min, niay be estimated simply by comparing the bremsstrahlung cooling rate to the total heating 
rate; we find 

^h,min — CoTp O (31) 

where Co ^ 0.031 and q = 1/2 when electron-proton bremsstrahlung dominates, and Co ~ 0.16 and 
q = 1 when electron-electron bremsstrahlung dominates. These estimates are confirmed by the more 
careful calculations shown in Figs. 6, 7, and 8. 

Where Ih < lh,min, we expect © oc ij/''. The curves in Figs. 3, 4 and 5 allow us to qualitatively 

check this behavior. Because q is itself an increasing function of O for the interesting range of 
temperatures, we expect Q to increase most steeply with Ih in the bremsstrahlung-domiriated regime 
when is least. That occurs for the larger values of Tp and h/lh, an expectation which is at least 
qualitatively borne out in these figures. 

We also observe, comparing the cooling rate plots with the corresponding temperature and 
optical depth equilibrium figures, that I a becomes greater than Ib at very close to the same place at 
which z becomes > 1. 
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4-4-^- Thermal equilibrium in the inverse Compton regime: the luminosity 

enhancement factor 



We have now identified the range in hard compactness corresponding to conditions in which 
thermal Comptonization controls the cooling of the system. In this case, we have 

Ihc = Ih, 

and the considerations discussed in §2.2 should in principle apply. Except at the very highest 
Ih, wherever inverse Compton cooling dominates, externally-supplied photons dominate the 
internally-created. In this section we will discuss how one might attempt to analytically define scaling 
relations between Is, Ih, Q, and tt in this regime. 

When inverse Compton on externally-supplied photons dominates, we can evaluate the hard 
luminosity, Lh, as 

Lh = "I'-eC / n^^^x)xdx — Ls{l — e"'^'^), (32) 



3 tesc J 

where the second term is subtracted in order to obtain the net Compton luminosity, and nsc{x) is 
given by eq. (15). We define 

<x>,.X^4^, (33) 
J n8c{x)dx 

where the sufflx 1 means that, with this definition (see eqs. (15) and (20)), the average over photon 
energy includes only those photons which scattered at least once. These definitions yield a simple 
relation between the input and output compactnesses, which is valid when the spectrum of injected 
photons is much narrower than the output spectrum (an assumption well-satisfied here): 



Ih/ls 



3il+g^TT/3) 



(l-e-^n, (34) 



where (x/xo)i is the mean factor by which the energy of incoming photons is amplified before escape, 
the factor 1 -|- grrr/S is the ratio between the escape time for the photons and the time R/c, and the 
factor of e~'^^ corrects for those soft photons which pass through without scattering at all. 

As in Dermer et al. (1991), the compactness ratio can be formally expressed in terms of the 
luminosity enhancement factor, r;: 



r 47ri?3 2 

J^hc = 7: TUeC 



j x'hsoft{x')\ri{x' ,Q,tt) -I], (35) 



where hsoft{x') is the spectral soft photon injection rate, and rj is in general a function of 
the initial soft photon energy x' , in addition to its dependence on Q and tt- In our case 
nsoft{x') = no5(x' — xo)/(R/c). The integrand is multiplied by r? — 1 rather than r; so that Lhc is the 
net Compton cooling. Consequently, for the present problem, from eq. (35), we have 

7^ = 7^ = '7(a;o,6,rT) - 1, (36) 

and, comparing with eq. (34), 

^ = e-^ + | j^^V (1-e-^). (37) 

One commonly followed approach to computing r] is to adopt a simple model of inverse Compton 
scattering in which to reach a given amplification a/ao always requires exactly n{x/xo) scatterings. If 
X <^Q, n{x/xo) — ln(a;/a;o)/ln(^r) (Ar is defined in equation (7)). In terms of this model, 

oo 

(-> = VP„-(n), (38) 

Xo ^-^ 
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where Pn is the probabihty of escape after exactly ri scatterings. Note that the relation between 
(x/xo), accounting for those photons that are not scattered at all, and the average we have introduced 
above (eq. (33)), (x/xo)i, is simply 

{±) = {±),(l-e-^^) + e-^^. (39) 

Xo Xo 

Several different approximation schemes have been proposed to estimate this sum. The simplest 
{e.g. Zdziarski et al. 1990) effectively replaces Pn with S„m, where Sij is the Kronecker delta, and m 
is the mean number of scatterings before escape. Loeb et al. (1991) suggested approximating P„ by 
the Poisson distribution. For the problem at hand, however, P„ is best described by the solution of 
the "gambler's ruin" problem (see, e.g. Feller 1967). The high-n tail of this distribution falls rather 
more slowly than in the Poisson distribution. At the same time, the amplification A" increases 
exponentially with n. Therefore, the largest contribution to the sum in equation (38) is likely to 
come from a range of n rather larger than m, and the associated amplification can be larger by a 
very substantial factor. To check just how bad an approximation the Poisson distribution is, we have 
compared the predictions of the Loeb et al. model (using j/r as given by eq. (9)) to the actual values 
of r){Q,TT) found in our solutions. The ratio between the two exhibits a very large scatter, with a 
general offset in the sense that the Poisson model substantially underestimates the true amplification. 

For these reasons, an alternative description of Comptonization was proposed by Dermer et al. 

(1991). They basically suggested that the amplification factor (x/xq) should be computed directly 
from the output spectrum (which they evaluated with a Monte Carlo simulation). We approximate 
the output spectrum analytically by a power- law plus a Wien component, so we write 

(±),=f^M+fw^ (40) 

Xo Xo Xo 

where fpi , fw are the fractions of the scattered photons going into the power- law and the Wien 
component, respectively, and fpi + fw ~ 1- The power-law fraction is 

J°° npi{x)dx I 

where we have used equations (19) and (26). The average photon energies are 

<^.^>=h<P^"^^^'-''^ forO<a<l^ 
V^oj^ fora>l 

(here we neglect to specify the case for a = 1 exactly); and {xw) = 30. 

In fact, over almost the entire range of h/lh we have studied, fpi ~ 1, i.e., /„, << fpi. f^ 
increases with Ihc/h, but the greatest value we find for it is ^ 4 x lO^'^ for Ihc/h = 10^. At this 
extreme value, the energy in the Wien term begins to be competitive with the energy in the power-law 
term; everywhere else it is negligible. Thus, the power-law portion of the inverse Compton losses is 
nearly always the dominant one. 

We therefore make the approximations that fpi ~ 1 and fpi{xpi) ^ fw{xw). In this limit, when 
< a < 1 (nearly all our equilibria fall in this range). 



^ ~ (1 - e-"^) 



(4/3)ar(l -a)(e/xo)i- 



(43) 



Because a depends in a rather complicated way on Q and tt, there appears to be no simple analytic 
relation describing the thermal balance in this regime. 

On the other hand, when a > 1, the analogous relation becomes 

"(4/3)a/(a-l) 



1 + 5ttTt/3 



Once again, the absence of a simple functional relation between a and and tt prevents further 
reduction of this form. 
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4.4-3. Empirical scaling 



Our numerical results do, however, reveal several simple scaling laws. Fig. 9 shows O as a 
function of the corresponding equilibrium tt, for all our equilibria for which inverse Compton cooling 
dominates. It is apparent that the points for each value of h/lh describe very nearly a straight line 
in this diagram, and these lines are all very nearly parallel, with logarithmic slope ~ — 1. There is no 
break apparent in any of these lines, at tt = 1 or anywhere else. For fixed Tp, the spread in tt is 
caused by increasing numbers of pairs. At the same time, differing values of Tp which share the same 
la/lh and happen to yield the same tt also coincide in Q, demonstrating that tt is the appropriate 
physical variable for this correlation. 

We are thus led to plot, in Fig. 10a, the product (6tt) as a function of the ratio h/lhc, for 
the same choice of equilibrium solutions; Fig. 10b shows the corresponding values for a, again versus 
Is/lhc- Again, the relationships are very nearly linear, and with very little scatter. Fitting to the 
points in Fig. 10a yields the relation 

he/ Is (lOerr)'', (44) 

where /3 = 4 ± 0.5. This relationship is a much better description of Comptonization equilibrium than 
any expression of the form e" . This empirical scaling law is, of course, what lies behind the trends we 
described in §4.3. 

The physical origin of the correlation with Is/lh can be easily seen, at least qualitatively. The 
cooling rate depends on the seed photon energy density, which rises with increasing Is/lh- However, 
because even unsaturated Comptonization can significantly increase the photon energy density above 
the energy density directly injected, there is a nonlinear dependence on the two quantities which 
promote Comptonization, namely Q and tt- That the nonlinear dependence should be describable so 
simply is surprising, given the arguments of the preceding section, but it does appear to be correct. 

This scaling law is not completely new. Haardt & Maraschi (1991) demonstrated that h/lh fixes 
a function of 6 and tt, and in Haaxdt & Maraschi (1993) they presented results showing that when 
TT < i, this function depends approximately on the product Qtt rather than on the two variables 
separately (at least for O not too large). What is new here is two things: the extension of this relation 
to Tt > 1, and the scaling with h/lh- 

In addition, we find that the coefficient in this scaling law depends weakly on geometry. For fixed 
Is and Tt, a plane-parallel slab has a higher average soft photon energy density than a sphere because 
the mean number of scatterings to escape is greater. Greater photon energy density leads to stronger 
cooling, and hence, for fixed Ih, a lower equilibrium temperature. Consequently, the coefficient of Qtt 
in equation (44) changes from 10 to ~ 16. 

We have already remarked (§4.3) that in the inverse Compton dominated regime, a is essentially 
a function of h/lh alone. In Fig. 10b the plot shows that functional relation. Fitting to these points 
yields 

(fe) ■ 

Comparing the scaling laws shown in the two parts of Fig. 10, we find the very simple relation 



1.6 1 



ac.—, (46) 

as can be visually inferred from Figs. 10a, and 10b. The coefficient 0.16 applies in the case of spherical 
geometry; in slab geometry, it becomes 0.1. 

When Is/lh <C 1, there is hardly any scatter around the relation given in equation 45; however, 
when Is/lh is more than a few tenths and tt < 1, there is also a weak secondary dependence on tt of 
the form a ~ I.Att^'^ ■ For smaller h/lh or larger tt, the spectral index is essentially independent of 
Tt- 
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In slab geometry, the increased number of scatterings for fixed tt compensates for the lower 
temperature, so that the relation between spectral index and h/lh is essentially identical to the 
spherical case. Slightly more energy is carried off in the Wien component however, to make up for the 
lower energy cut-ofF in the power-law component. 

The scaling shown in equation 46 is in fact a fairly good approximation to the expression for 
a{Q,TT) given in §3.1.1 in the case of modest optical depth (1 < tt < 3) and low temperature 
(e <, 0.3): 

-ln[l-3/(4TT)] 3 

" 49 ^ 16(rTe) • 

This expansion predicts the correct scaling and very nearly the correct coefficient (0.19 instead of 
0.16). That the scaling law is in fact more broadly applicable is due to the fact that we consider only 
the portion of the Q - t plane in which equilibria having h/lh in the range 10~^ - 1 are found. This 
condition restricts Qt to be ^ 0.1, and when this condition is met, equation (46) continues to be 
a fairly good approximation to the exact expression even when r < 1. That is, equation (46) only 
applies when the plasma is in thermal balance. 

,]ust as for the temperature scaling law expressed in equation (44), hints of this relation between 
spectral index and h/lh have also appeared in the literature. Gil'fanov et al. (1994) argued on a 
qualitative basis that a should depend primarily on h/lh, and vary in the sense that we have found. 
Working solely on the special case h/lh — 0.5, Haardt & Maraschi (1993) computed the spectral index 
over the range 0.01 < tt < 1 and found that it hardly varied (the sharp gradient in a they found as 
Tt approaches 1 is due to a breakdown in their approximations in this regime) . What we have done is 
to confirm that this relation exists, and to discover its quantitative character. 

Wc have now arrived at a much more powerful description of the scalings between physical 
parameters in the inverse Compton-dominated regime. It is no accident that they have been found 
by a systmatic approach focusing on causal parameters, rather than fitting to phenomenological 
parameters. In fact, an approach which treated Q and tt as independent parameters would have 
been incapable of discovering these correlations because there would have no sensible way to limit the 
region of the © - tt plane considered. Only framing the problem with Ih and h/lh as the independent 
variables reveals these correlations. 



5. Consequences for Observable Systems 



5.1. The Natural Ranges for 1^, Is/lm ^iid Tp 



5.1.1. Ih 



As we have just seen, the temperature is confined to a rather narrow range over a very broad 
span of Ih and h/lh, provided 0.1 ^ Tp ^ /ew. However, this volume of parameter space, large as it 
is, is not infinite. We now argue that it would be however, very unlikely for any accreting black hole 
to be found outside this volume. 

First, as many have previously pointed out {e.g. Lightman and Zdziarski 1987), if the 
X-ray-emitting plasma is roughly spherical, Ih may be rewritten in terms of quantities more directly 
related to the accretion process: 

Lh L 27r^ie/me , 

i-h = ^j ] , (47) 

L Le r Tg 
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where L/Le is the total luminosity relative to the Eddington luminosity, fXe — l-2mp is the mass per 
electron, and r/rg is the size of the region relative to the gravitational radius of the central mass. 
Most of the gravitational energy release in accretion onto a black hole takes place around 5 - lOvg, so 
it would be difficult to make Ih more than 3 x 10^. On the other hand, if the portion of the dissipation 
going into the X-ray-emitting plasma in that region is at least ~ 10~^Le, Ih is > 10~^, the lower limit 
of the region in which we find almost constant temperature. Thus, if there is enough energy going 
into the X-ray-emitting plasma near a black hole to make it detectable, Ih will fall into the constant 
temperature range. 

In this respect, coronae near the accretion disks around weakly-magnetized neutron stars should 
be very similar to those near black holes. Accreting white dwarfs genorically should have Ih at least 
~ 10~^ times smaller for equal mass accretion rate, and are therefore likely to fall just below the 
constant temperature range of Ih even if Lh/L ~ 1. With Ih ~ 10"'' 10"'\ corona; in the inner 
regions of accretion disks around white dwarfs are likely to be in a bremsstrahlung-dominated regime 
(see Fig. 6) with temperatures of a few keV. 

The assumption of roughly spherical geometry is not entirely innocuous. It is almost certainly 

valid for largo Ih, where the plasma is always pair-dominated. When that is true, the electron thermal 
speeds are comparable to the escape speed even for the smallest possible r/rg, and the positrons 
cancel out any restraining electric fields. In addition, the effective Eddington luminosity is reduced 
by - 10-\ How ever, at the small Ih end, whore the plasma is not pair-dominated, a less symmetrical 
geometry {e.g. disk-like) can make a difference. For fixed tt, the ratio of bremsstrahlung luminosity 
to inverse Compton luminosity is oc {r/h)^lf^^, where h is the smaller dimension. Thus, flattened 
geometry can cause the lower bound of the constant temperature range of Ih to increase substantially. 



5.1.2. h/lh 

If Ls^intr/Lh -C C, Is/lh is Controlled by essentially geometrical factors and is (almost) always 
< 1 (see equation 29) . A common assumption is to suppose that the hot plasma forms a slab bounded 
on one side by a (comparatively) cool accretion disk (Liang 1979, Haardt & Maraschi 1991). In that 
case, Is/lh — 0.5, with some order unity corrections due to anisotropy in the upscattering (Ghisellini 
et al. 1991; Haardt 1993). When the intrinsic soft luminosity is small, while it is very difiicult to make 
Is/lh > 1, there are many ways it could be <^ 1. Both C and r^/d^, for example, must be less than 
unity, and could be much less. In addition, modest relativistic motion can alter the effective value of 
Is/lh because this ratio should be evaluated in the bulk frame of the hot plasma. If the hot plasma 
moves away from the soft photon source with speed /?, the comoving energy density of soft photons 
falls as [7(1 -|- /3)]~*. Motion towards the soft photon source would have the opposite effect, of course, 
but it would take rather contrived dynamics to produce such a situation. 

When Ls,intr is not negligible, h/lh can easily be > 1. However, as shown in §4.4.3, the 
equilibrium temperature on the flat part of the curve is oc {ls/lh)~^^'^, so rather large values of h/lh 
are necessary in order to substantially cool the plasma. 



5.1.3. Tp 

As is demonstrated in Figures 3, 4, and 5, the total optical depth tt becomes independent of Tp 
(and a function only of Ih and h/lh) when Ih is large enough for the number of pairs to be much larger 
than the net lepton number. Regulation of Tp is therefore only relevant at smaller Ih- The remainder 
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of this subsection therefore focusses on the range of parameters in which there are relatively few pairs 

and tt — Tp. 

So far we have relied ou a mental picture in which the hot X-ray-emitting plasma is physically 
well separated from the cool source of soft photons. If they actually touch each other (which is a 
statement almost equivalent to h/lh ~ 1), we must consider thermal mixing via electron conduction. 
Good magnetic connectivity is probably a reasonable assumption because the most likely way to 
actually supply the energy to dissipate in the hot plasma is through the passage of MHD waves from 
an accretion disk to an adjacent hot corona (Stella and Rosner 1984; Tout and Pringle 1992). If 
wave-particle scattering does not grossly diminish the hot electron mean free path, it makes little 
sense in this case to speak of a hot plasma unless it is significantly thicker than a Coulomb mean free 
path. If it were otherwise, much of the volume occupied by the hot electrons would actually be shared 
with the cool matter (which would then be evaporating into the hot gas: Balbus and McKee 1982). 
This restriction places a lower limit on tt, 

rT>^eV In A, (48) 

o 

where In A is the usual Coulomb logarithm, and b changes from 2 to 4 as increases from -C 1 to 
^ 1. If O > 0.1, the lower limit on tt due to thermal conduction is ~ 0.01. 

When the optical depth becomes large, gradients in the radiation pressure develop. Large forces 
can be associated with these gradients; the minimum radiation force is larger than the force due to gas 
pressure gradients when Ih > AtvQ in spherical geometry. In the presence of gravity, radiation pressure 
dominance often leads to dynamical instabilities {e.g. Arons, Klein, & Lea 1987; Krolik 1979). We 
speculate that these dynamical instabilities may place an upper limit on the optical depth by causing 
the hot plasma to break up into smaller pieces whose individual optical depths are no bigger than a 
few. 



5.2. The Connection between Spectral Index and Intrinsic Soft 
Luminosity or Source Geometry 

AGN exhibit a rather small range of power law indices in the X-ray band. Type 1 Seyfert 
galaxies AGN have intrinsic spectral slopes in the 2-20 keV band of ~ 0.9 with very little scatter 
if the effects of Compton reflection and partially-ionized absorption are factored in (Mushotzky 
et al. 1993). Higher luminosity AGN (i.e. quasars) have similar spectral slopes, although possible 
contributions of reflection and ionized absorbers are harder to determine (Mushotzky et al. 1993). 

If the thermal Comptonization model is correct, we can infer that h/lh is most likely held stably 
at ~ 0.1 in all of these AGN. Such a small value of Is/lh implies that at least one of the following 
two conditions applies: both Ls,intr/Lh and the covering fraction of any reprocessing surface must 
be small; or the distance from the source of soft photons must be rather greater than the size of the 
X-ray emitting region. In type 1 Seyfert galaxies, the frequent evidence for strong Compton reflection 
suggests C ~ 0.5; where this is the case, {R/d)^(f> must be ~ 0.1. Note that contributions to C from 
the obscuring torus (Krolik et al. 1994; Ghisellini et al. 1994) automatically have {R/c[f -C 1. 

Moreover, because a may be used to infer la/lh, and la/lh fixes Qtt, a characteristic a ~ 0.9 
implies that Qtt ^ 0.1 - 0.2. The smaller value obtains when the plasma is a slab, the larger value 

when it is spherical. 

The only possible alternative to these inferences is a very specific version of the coronal model 
{e.g. Haardt & Maraschi 1991, 1993; Zdziarski et al. 1994), in which C ~ 0.5 and J? ~ d. Because 
Is/lh — 0.5 in these models, they fall in the regime in which there is a small secondary dependence 
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of a on tt- For this value of h/lh, the observed spectral slope can be produced if tt is somewhat 
less than 0.1. The uniformity of AGN spectra would then require a rather narrow distribution of 
Tt- It follows as a corollary that would also have a narrow distribution, centered on ~ 0.8. Heat 
conduction (§5.1.3) might pose a physical self-consistency problem for models of this sort. 

If h/lh is really ~ 0.1, because the observed ratio of Ls/L^ is often order unity or somewhat 
greater, some moans must be found to reduce the local value of h/lh in the the X-ray production 
region. For example, it might not be immediately adjacent to the region where the bulk of the 
ultraviolet luminosity is radiated. Another possibility is that the X-ray region moves away from the 
disk fast enough that the seed photon intensity is diminished by redshifting. All these schemes are 
made easier if the intrinsic soft luminosity is minimized. 

In stellar black hole binaries, the X-ray spectrum is most simply described as the sum of a 
quasi-thermal soft component and a hard power-law (Tanaka 1989, White 1994). The soft (i.e. T ~ 1 
- 3 keV) component is presumably due to dissipation in the accretion disk proper, where the density 
and optical depth are great enough to approximately thermalize the spectrum. Weakly-magnetized 
accreting neutron stars should also possess such a soft component, although much of the dissipation 
may occur in shocks near the surface of the neutron star rather than as a result of gradual viscous 
dissipation in the disk proper. In addition, accreting neutron stars may have rather weaker hard 
components than accreting black holes if the disk is interrupted more than a few neutron star radii 
away from the surface. 

The hard power-law seen in accreting black hole systems is generally thought to be the result 
of thermal Comptonization (Shapiro et al. 1976). Spectral indices of this hard component exhibit a 
greater range than in AGN; examples are known from a ~ 0.3 to a ~ 1.5 (Tanaka 1989; Ballet et al. 
1994; Gil'fanov et al. 1994), but the center of the distribution is fairly close to the average spectral 
index in AGN. We therefore expect a larger range in Is/lh, but the mean inferred h/lh for stellar black 
holes is also ~ 0.1. As for AGN, very optically thin coronal models can also produce spectra with 
a ~ 1, but it may be difficult for them to generate spectra as hard as a ~ 0.3. 



5.3. Coupling between Fluctuations in the Luminosity and the Spectrum 

The simplest model for X-ray fluctuations one might imagine is one in which the geometry and 
the relative division of heating between hot plasma and cool gas {i.e. the factors which determine 
Is/lh) are fixed. If the column density of the corona is fixed independently, changes in luminosity 
drive changes in Ih, but not in the other parameters. In the large volume of parameter space in which 
the temperature varies only slowly with Ih, one would then expect that luminosity fluctuations should 
raise the general continuum level without altering its spectrum, for © and tt, which determine the 
X-ray spectrum, are nearly invariant. 

In only one AGN (NGC 4151) are the effects of absorption and Compton reflection weak enough 
that one can attempt to measure possible variations in the intrinsic continuum shape (Yaqoob 1992). 
In this object (whose mean spectral index is comparatively small, ~ 0.5), a rises from 0.3 to 0.7 while 
the 2-10 keV flux grows by about a factor of four. In the framework of our model, this behavior 
would be interpreted as an increase in la/lh from ~ 0.001 to ~ 0.04. Whether this change in Is/ln 
can be genuinely associated with a change in total luminosity is hard to say because the total power 
in spectra this hard is dominated by the highest energy photons, and the data which exhibit this 
correlation say nothing about the flux at the photon energies which dominate the output. 

There is a much greater wealth of observations regarding stellar black hole X-ray variability 
than regarding AGN X-ray variability. In the case of stellar black holes, the amplitude of the hard 
component often fluctuates by very large amounts: changes of three orders of magnitude have been 
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seen. Yet despite these very large fluctuations, in most examples the spectral index has remained 
essentially fixed (Tanaka 1989; Churazov et al. 1993; Gil'fanov et al, 1994). If the ratio of soft flux to 
hard were constant, this behaviour would, of course, be precisely in keeping with the simplest model 
one could imagine. 

However, the hard and soft components in fact appear to vary quite independently. Therefore, if 
the power law index of the hard component remains constant despite changes in L,,/Lft,. there must 
be an additional element in the model beyond the simplest case. One possibility is that the main 
source of soft photons is quite distant from the radiation region for the hard component, so that 
the soft photons serving as Comptonization seeds are entirely from local reprocessing with constant 
reprocessing fraction. Another possibility is that the distance between the Comptonizing plasma and 
the main source of soft photons changes in such a way as to cancel the change in luminosity ratio. 
For example, if the Comptonizing plasma is heated by dissipation of MHD waves originating in an 
accretion disk, an increase in the intrinsic soft luminosity of the disk might be accompanied by an 
increase in the flux of MHD waves into the hotter plasma. As a result, the hot plasma could be pushed 
farther away from the disk, so that the soft photon density in the Comptonizing region is reduced. 

In at least one example (Nova Muscae: Ebisawa et al. 1994), the spectral index of the hard 
component did change. In this the soft component faded, the spectral index of the hard 

component fell from ~ 1.5 to ~ 0.7. At least in a qualitative sense, this follows the general trend of 
decreasing la/lh leading to harder Comptonized power-laws. 

The direct relation between spectral shape and h/lh breaks down on the shortest timescales. 
Spectral equilibration requires a time teq ~ ln{Q /xo)R/ {ctt) because the highest energy photons are 
scattered ~ ln(0/a;o) times. On the other hand, Is/lh could easily change in a time J?/c, which is 
typically an order of magnitude less than teq- Therefore, these correlations should not be followed on 
timescales shorter than teg- 



6. Summary 



In this paper we have demonstrated that the prevalence of ~ 100 keV thermal Comptonization 
spectra in the X-ray emission from accreting black holes is no coincidence. Over a very wide range of 
dimensionless heating rate (i.e. compactness Ih ), and normalized strength of seed photon supply (i.e. 
Is/lh), and a somewhat narrower range of net lepton optical depth (rp), thermal Compton scattering 
dominates the heat balance, and the temperature is nearly always approximately of this magnitude. 
In fact, for fixed Is/lh and Tp, the temperature is virtually independent of heating rate over a dynamic 
range in Ih that is generally at least several orders of magnitude. 

After quantifying the boundaries of this constant temperature regime, we have also argued that 
nearly all accreting black hole systems should fall within its (large) volume of parameter space. That 
Ih and Is/lh should fall within that range seems very likely, given what we know of black hole accretion 
dynamics; that Tp should be in the range giving this result is less certain, but we have raised several 
plausible arguments toward this end. The compactness is determined by the nature of accretion into a 
relativistically deep potential; the seed photon supply follows simply from geometrical considerations; 
the optical depth may be determined by a combination of thermal conduction and radiation pressure 
dynamics. 

We have also shown that simple scaling laws relate the temperature, total {i.e. including pairs) 
Thomson optical depth, and the seed photon supply h/lh- h/ln (almost) uniquely determines the 
product Gtt, and, to roughly the same level of approximation, Is/lh alone determines the spectral 
index a of the power-law segment of the Comptonized spectrum. 
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On the basis of this relationship between a and Is/lh we have argued that Is/lh in AGN is 
most likely to be generically ~ 0.1. Such a small value of Is/lh requires any thermalizing surface (for 
example, an accretion disk) to be comparatively distant from the X-ray emission region, and suggests 
that the intrinsic (as opposed to reprocessed) soft photon luminosity is comparatively small. If instead 
the thermalizing surface is nearby, so that h/lh = 0.5 or more, the optical depth of the corona must 
be rather small, and its value is quite tightly constrained; the coronal temperature is then similarly 
tightly constrained. In stellar black hole binaries, the observed range of power law indices is larger 
than in AGN, so a larger range of h/lh is permitted, but this range is also roughly centered on 0.1. 
When stellar black hole binaries make transitions between "soft" and "hard" states, the slope of 
the Comptonized power law usually remains constant, even while the relative amplitude of the soft 
thermal and hard Comptonized components change by orders of magnitude. This fact requires the 
luminosity in seed photons fed into the Comptonizing plasma to not be simply proportional to the 
soft luminosity we observe. 

Wc thank Francesco Haardt and Andrzcj Zdziarski for many enlightening conversations. Wc also 
thank Andrzcj Zdziarski for much guidance in the nuances of pair equilibria and Comptonization. 
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Figure Captions 

Figure 1 A global view of our equilibria, seen as functions of Q for Tp = 1.0 and xq = 10~^. 
The various curves are each for a different fixed value of no /rip (and, correspondingly, of Is', solid is 
no/Up = 0; dotted is no/rip = 100 (i.e., h = 3.14 x 10"^); short dash is no/wp = 10^ (h = 3.14 x 10"^); 
long-dash is no/rip = 10*(Zs = 0.314); dot-short-dash is no/rip = 10^ {Is = 3.14); dot-long-dash is 
no /rip = 10® (Is = 31.4). In the first three panels Ih increases upwards; in the fourth it is the opposite. 

a) The normalized pair density n+/np. Pairs dominate only when the equilibrium is well on the 
high-lh side of the temperature maximum, but they can significantly influence the temperature even 
when z ~ 0.1. 

b) The total Thomson optical depth tt- 

c) The heating rate in terms of compactness Ih- The curves for different no /rip converge when 
they become pair-dominated. 

d) The power law index of the Comptonized output spectrum. No curve is plotted for no/rip = 
because this spectral component does not exist in that limit. 

Figure 2 The relative importance of different cooling processes as functions of Ih for the specific 
case Tp = 1, no/rip — 10* {h = 0.3). The solid line is inverse Compton cooling on external photons; 
the dotted curve is inverse Compton cooling on internally created photons; the dashed curve is 
annihilation cooling; and the dot-dashed curve is bremsstrahlung. For most of the parameter range, 
inverse Compton cooling on externally-created photons accounts for nearly all the cooling, but at the 
highest values of Ih, it is overtaken by inverse Compton cooling on internal photons. 

Figure 3abcd The principal physical quantities as functions of Ih for Is/lh = In this figure and 
Figs. 4 and 5, we show O (panel a), tt (panel b), a (panel c), and yr (panel d) as functions of Ih- In 
each of Figs. 3, 4, and 5, the solid line corresponds to Tp = 0.1; the dotted line to Tp = 0.5; the short 
dashed line to Tp = 1; the long dashed line to Tp = 2; and the dot-dashed line to Tp = 3. 

Figure 4abcd The principal physical quantities as functions of Ih for Is/lh ~ 0.1. The different 
curves correspond to different values of Tp, as described in the caption to Figure 3. 

Figure 5abcd The principal physical quantities as functions of Ih for Is/lh = 0.01. The different 
curves correspond to different values of Tp, as described in the caption to Figure 3. 

Figure 6abcd The relative contributions of the different cooling mechanisms for l^/lh = 1- The 
solid line is bremsstrahlung (the sum of e-e, e-p and e-e"*"); the broken line is total inverse Compton 
(on both internally and externally created photons); the dotted line is annihilation cooling. In some 
cases annihilation cooling is so weak that it accounts for less than 10~® of the total, and its curve does 
not appear. 

a) Tp = 0.1 

b) Tp = 0.5 

C) Tp = 1 

d) Tp = 2 

Figure 7abcd The relative contributions of the different cooling mechanisms for Is/lh = 0.1. The 
mechanisms associated with each curve are as described in the caption to Figures 6. 

a) Tp = 0.1 



b) Tp = 0.5 
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C) Tp = 1 

d) Tp = 3 

Figure 8abcd The relative contributions of the different cooling mechanisms for Is/lh = 0.01. The 
mechanisms associated with each curve are as described in the caption to Figures 6. 

a) Tp = 0.1 

b) Tp = 0.5 

c) Tp = 1 

d) Tp = 3 

Figure 9 Temperature in the Compton-dominated regime as a function of tt for a variety of seed 
photon supply rates. The symbol coding is: Is/lh ~ 1 represented by filled circles; Is/lh ~ 0.5 by solid 
squares; Is/lh = 0.1 by crosses; h/lh = 0.01 by open circles; and h/lh = 0.001 by open triangles; the 
open square represents an equilibrium corresponding to h/lh = 10. 

Figure 10 The scaling of Comptonization model parameters with Is/lh- In these figures we have 
plotted results representative of all the models we computed. 

a) Qtt as a function of Is/lh- Despite differences in Ih and Tp, all models fall very close to a 
straight line in log(0rT) vs. \ogls/lh- For fixed Tp, the spread in tt is due to pair creation. 

b) a as a function of h/lh- Although three independent parameters are required to specify these 
models {Ih, h, and Tp), the spectral index of the power-law component of the output spectrum depends 
almost solely on one parameter, h/lh- When h/lh — 1, there is a weak secondary dependence on tt- 



